Analysis of f-p model for octupole ordering in NpC>2 
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In order to examine the origin of octupole ordering in Np02, we propose a microscopic model 
constituted of neptunium 5/ and oxygen 2p orbitals. To study multipole ordering, we derive effective 
multipole interactions from the f-p model by using the fourth-order perturbation theory in terms 
of p-f hopping integrals. Analyzing the effective model numerically, we find a tendency toward 
antiferro-octupole ordering. 
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In the research field of condensed matter physics, it 
is one of currently important issues to unveil the mech- 
anism of multipole ordering phenomena frequently ob- 
served in /-electron systems. It has been widely rec- 
ognized that quadrupole ordering realizes in several /- 
electron compounds, but recently, in addition to dipole 
and quadrupole ordering, a possibility of ordering of 
higher order magnetic multipoles, i.e., octupoles, has 
been discussed intensively. 

A typical candidate with octupole ordering is the low- 
temperature ordered phase of NpC^, 1 since time reversal 
symmetry is broken in this phase^^ but the detected 
internal field is too weak to be ascribed to dipole order- 
m g ( 2 J 3 i 4 j n( jeed, several experimental facts^i^ can be 
consistently explained by assuming longitudinal triple- 
q octupole ordering; 9 In addition, recent experi- 
ments on the 17 O NMR also support the triple-q ordered 
state3ii 

In order to understand why such higher-order multi- 
pole order is realized in NpC>2 , it has been highly required 
to proceed to the microscopic research. In general, it 
is difficult to develop a microscopic theory for complex 
multipole ordering in /-electron systems beyond the phe- 
nomenological level, but it has been recently proposed to 
construct a microscopic /-electron model on the basis of 
a j-j coupling scheme. 12 Following this proposal, we have 
studied the /-electron model on an fee lattice composed 
of Np ions with hopping integrals via (ffcr) bonding, and 
actually found the triple-q octupole ordering^ However, 
as shown in Fig. [fl^a) i oxygen anions exist between Np 
ions in actual crystal structure. Thus, it is important to 
clarify how the octupole ordering in the /-electron model 
is affected by oxygen anions. 

In this paper, we construct a more realistic model in- 
cluding also p orbitals of oxygen anions in addition to / 
orbitals of Np ions. We derive an effective multipole in- 
teraction model by evaluating the exchange of / electrons 
via p orbitals within the fourth-order perturbation with 
respect to p-f hopping integrals. By analyzing the effec- 
tive model, we again find a tendency toward the triple-q 
octupole ordered phase, indicating that the /-electron 
model on the fee lattice have grasped the essential point 
on the appearance of octupole ordering of NpC>2 . 

Let us discuss local /-electron states of actinide diox- 
ides based on the j-j coupling scheme, in which we first 



FIG. 1: (a) Crystal structure of Np02- a is the lattice con- 
stant, (b) 8-site cluster (gray spheres) on the fee lattice taken 
in our numerical calculation. 



include the spin-orbit interaction and consider only the 
states with total angular momentum j—5/2. From a 
quantitative viewpoint, this simplification may not be ap- 
propriate for actinide dioxides, since they are insulators, 
and the LS coupling scheme is expected to work well to 
describe /-electron states in these materials. However, 
both schemes are continuously connected to each other 
by changing the ratio of the strength of the spin-orbit 
interaction and the Coulomb interaction, as long as the 
symmetry of the ground state is not changed^ Thus, 
on the basis of a spirit of adiabatic continuation, we ex- 
pect that qualitative properties at low temperatures can 
be grasped whether we choose the j-j coupling or the 
LS coupling schemes as a starting approximation. The 
,7=5/2 states split into Tj doublet and T% quartet un- 
der a cubic crystalline electric field (CEF). Since the Tj 
wavefunction extends along the [111] direction and an 
oxygen anion locates in this direction, the level is ex- 
pected to be higher than the Tg level. If we assume that 
the level splitting A between Tg and Tj is large enough, 
CEF ground states for f 2 , f 3 , and / 4 are obtained by 
accommodating two, three, and four electrons into Tg 
levels, leading to T§, Tg, and Fi, respectively, consis- 
tent with experimental results for UO21 15 Np02,^ and 
Pu02^ respectively. Thus, we ignore the T 7 state to 
discuss the ground state of actinide dioxides in the j-j 
coupling scheme. 

We have three comments on the CEF level scheme, (i) 
In our picture, the first excited state of Pu02 should in- 
clude three Tg and one T 7 electrons, indicating that the 
excitation energy provides the lower limit for A. Since 



2 



the CEF excitation energy of PuC>2 is 123 meVji£ A 
should be larger than 1400K, consistent with the initial 
assumption, (ii) The / 3 state in NpC>2 is regarded as the 
one-hole state in T%. In the following, we use a hole pic- 
ture and "electron" denotes such a hole, (iii) Among the 
f 2 states, the Ts triplet is the ground state as observed 
in UO2. Since the CEF excitation energy in UO2 is as 
large as 150 meV^ we consider only the triplet among 
f 2 intermediate states to study exchange processes of / 3 
ions in Np02- 

The Tg quartet consists of two Kramers doublets, 
and it is convenient to introduce orbital index r (= 
a, (3) to label the two Kramers doublets and spin in- 
dex a (=1, i) to distinguish the two states in each 
Kramers doublet. In the second-quantized form, an- 
nihilation operators for electrons are defined by 

/rat =\/5/6 a r5/2 + W^ a v-3/2, /raj. =V^76 fl r-5/2 + 

W® a r3/2, / r/ 3T=a r i/ 2 , and /r/3i=a r -i/2, where a r]z is 

the annihilation operator for an / electron with the z- 

component j z of the total angular momentum j=5/2 at 

p 

site r. Multipole operators are usually expressed as X r T , 
where T 7 denotes symmetry. The explicit forms of X r '' 
in the Tg subspace are found in Ref. ll3L 

Now we show the f-p model for Np02, given by 
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H = Hf +Hp + Hkin, 



(1) 



where TLf and TL P are the local /- and p-electron terms, 
respectively, and TLun denotes the hybridization term 
between / and p electrons. Among them, the local /- 
electron term is explicitly given by 



Np 



(Np) 



4--.(d)*--2 



(c) 



(Np 



i*-4 2- 



Np 



(d) 



(Np. 



(e) 



(OU--2 



Np 



Np 



Np 



n r j3 



FIG. 2: Fourth-order processes with respect to p-f hoppings. 

(a) Electrons doubly occupy a Np site. Electrons occupy (b) 
two or (c) one O site simultaneously, and then exchange. Elec- 
trons occupy (d) two or (e) one O site simultaneously, and 
then return to the initial sites. The solid and dashed lines in 

(b) -(e) distinguish paths for different electrons. 



where v is a vector connecting nearest-neighbor Np and 
O sites and s symbolically denotes both spin and orbital 
of an / electron. The hopping integrals trier s are repre- 
sented in terms of (pfcr) and (pfir) by using the Slater- 
Koster tabled 

From the f-p model Q , we derive an effective interac- 
tion by using the fourth-order perturbation theory with 
respect to the p-f hopping integrals. Here we note that 
the fourth-order processes are classified into two types: 
(a) Two electrons doubly occupy a Np ion in the interme- 
diate states, as schematically shown in Fig.|2a). (b)-(e) 
Two electrons occupy one or two oxygen sites simultane- 
ously in the intermediate states, as shown in figures|2Ib)- 
(e) . The effective Hamiltonian is arranged in the form of 



J ^ ] fraafrficr' /raff' /r/3cr 
r,fT,ff' 

J' X! frrlfrTlfrr'lfrr'h 



(2) 



where e/ denotes the energy level of Tg, n TTa = f£ TtT f r T<j, 
and n TT =J2 a n rr<j- The coupling constants U, U' , J, 
and J' denote the intra-orbital Coulomb, inter-orbital 
Coulomb, exchange, and pair-hopping interactions, re- 
spectively. 

The local p-electron term is expressed as 



xm.cr 



(3) 



(on-site interactions for p orbitals), 



where e p denotes the energy level of p orbitals of oxygen 
ions and p xma is the annihilation operator for an electron 
with the z-component m of orbital angular momentum 
and real spin a at O site x. Finally, 7Ykin is given by 

Hkin = (Cff; S Pl+*mff/™ + h - c -)> (4) 
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(r,r'),si-S4 



^s 3 ,s 4 ;si,s 2 /rS3/rsi / r ' S4 /r's 2 > (5) 



where (r, r') denotes the pair of nearest-neighbor sites 
and the effective interaction / is given by 



(6) 



k— a— e 



In the following, we will briefly explain each term in- 
cluded in /. 

In order to derive multipole interaction from the 
process (a), we simply replace /-electron hopping i£ , 
in the effective interaction / of the /-electron model 13 
with the effective /-/ hopping T^ s , via p orbitals. For 
the purpose, we define /-electron hopping between sites 
r and r + pL via p orbitals of the oxygen ion at r + u, 
given by 

^s^s' ~ ( e f ~ e p) y ) fences) ^ma\s' • CO 
ma 

Then, we obtain T^ S ,=J2 U T^'Jf , since the effective /- 
electron hopping consists of two paths through different 
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oxygen sites. For instance, along the [110] direction, the 
effective /-/ hopping integral is given by j^/ 2 >°/ 2 > ) = 

T (a/2,a/2,0),(a/4,a/4,a/4) + r (o/2,a/2,0),(a/4,o/4,-o/4) whcrc 

a is the lattice constant. Using T^f , , we express 1^ as 

(■frir'-ru p* p in 1 
|V J s';s 4 -' r u-s 3 ,s' r u;si,s- L i 



r(a) /x 
S 3 ,S4;Si,S2 



u,s,5' 



nr — r 



(8) 



(rpT—T \* p* p 



rpr — T 



where P u - S ,s' denotes the inner product between one of 
the r$ triplet states denoted by u and the f 2 state labeled 
by s and s' . 

We note that the effective /-/ hopping integrals have 
the same form as those via (//cr) bonding^ for instance, 

T^ a J m « [trr'Sr* + C X a% T ,5 aa , + C^a^}, (9) 

where cr are Pauli matrices, and c\ and C2 are constants 
depending on (pfcr) and (pftv). This fact indicates that 
the form of the hopping integrals are restricted by /- 
electron symmetry, and the simple (//cr) model may 
grasp properties of actual materials with complex struc- 
tures. 

Concerning processes (b)-(e), we consider the effect 
of the Coulomb interaction at oxygen sites, symbolically 
expressed as ll Up" . In this paper, we study two limiting 
cases, U p =0 and U p =oo. For U p =0, multipole interaction 
contains all the processes (b)-(e), given by 
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J S3 ,S4',Sl ,52 
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J S3 ,S4;Sl ,52 




r(d) *x 
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= 2(e 



S 4 ;S1^S3;S2 ' 



and 



-1 \ * nv c-fJ.+l 

/ j S4;«i s 3 ;s 2 



j-(e) ti _ _r(d) A* 

i S3,s 4 ;si,s 2 S3,s 4 ;si,S2' 



respectively, where 



ma 



(10) 

(11) 

(12) 

(13) 
(14) 



Note that the sum of processes (d) and (e) merely 
becomes energy shift, since I, 



(c) It 

S 3 . s 4;Sl,S2 



r(d) A* 



1 S3,s 4 ;s 11 s 2 



<5 SlS3 (5 S2S4 , and such terms can be eliminated in the 
present discussion. After all, the multipole interaction 
for Up—0 is given by 



jfj, _ T(a) ti 

S3,s 4 ;si,s 2 s 3 ,s 4 ;si,S2 



-"«,. (15) 



Note that this term is expressed by only the effective /-/ 
hopping integral. 
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FIG. 3: Phase diagrams obtained from the multipole cor- 
relation functions. Here 3gF denotes Ts g moment with q = 
(0,0,0), and so on. The antiferro (AF) ordering vector is 
q = (0, 0, 1) in units of 2n/a for all the AF phases, (a) Phase 
diagram for U P =Q. (b) Phase diagram for U p =0 in a mag- 
nified scale around (U — J)/(e p — e/ + U — - J)=0.145 and 
(p/7r)=0. (c) Phase diagram for U p =oo. 



Now we move on to another limiting case U p —oo, in 
which the effect of U p is included by prohibiting processes 
(c) and (e). Thus, the effective interaction for U p —oo is 
given by 



S 3 ,S 4 ;S1 ,S2 



r(a) fj. 



r(b) ft 

S3,s 4 ;si,s 2 



. r(d) M 
J -s a ,s i ;s 1 ,s 2 ■ 



(16) 



It should be noted that in the process (b), a couple of 
electrons exchange their sites by avoiding the effect of 
Coulomb interactions. Such a term is characteristic to 
the crystal structure of Np02, and it has a tendency to 
stabilize the octupole ordering, as shown later. 
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Using the effective Hamiltonian, we evaluate nu- 
merically the multipole correlation function, ^q T — 

0-/N) J2 r .r> e lcl < r ~ r ">{Xr''X^), where (■ • • ) denotes the 
expectation value using the ground-state wavefunction. 
Here we take N=8, as shown in Fig. ^b). Figures E{a) 
and (c) show phase diagrams, presenting the multipole 
moment which has the largest value in the correlation 
function at each parameter set, for U p =0 and U p =oo, 
respectively. For U p —0, as shown in Fig. |3{b) in a mag- 
nified scale, there is a very small, but finite region of 
the antiferro-octupole [q = (0,0, 1) in units of 27r/a, 
5uAF] phase. For U p —oo, on the other hand, the re- 
gion of the 5uAF phase becomes much larger than that 
for U p =0. The 5uAF phase locates in the parameter re- 
gion with e p — €f <C U' — J, in which processes (b) and 
(d) work effectively. Since process (d) provides only a 
quadrupole interaction, we conclude that the stabiliza- 
tion of the 5uAF phase originates from the process (b). 

Note that the Coulomb energy in the f 2 intermediate 
state U' — J is expected to be in the order of 1 eV, but 
we cannot estimate the difference of the energy levels of p 
and / electrons e p — e/ within the present theory. We also 
note that in the phase diagrams, the 5uAF phase appears 
for (pfn) ~ 0. Since (pfa) and (pfir) are treated as pa- 
rameters in this paper, we have no clear answer why the 
absolute value of (pfn) should be so small for the appear- 
ance of the 5uAF phase. In order to confirm the reality 



of the parameter region for the octupole phase obtained 
in this study is realistic or not, it is highly requested to 
perform the band-structure calculations for NpC-2. This 
is one of future problems. 

In summary, on the basis of the f-p model for NpC>2, 
we have found a finite region of the 5uAF phase for both 
cases of U p =0 and U p =oo, Thus, we expect that this 
property is retained in the actual situation for Np02 with 
finite Up. While the ordered state cannot be entirely 
determined within the present small-cluster calculation, 
it is confirmed that among structures with q=(0,0, 1) 
and equivalent ones, the triple-q structure is energetically 
favorable^ since the T$ u moment in the T% subspace has 
the easy axis along the [111] direction. We emphasize 
that the T^ u antiferro-octupole phase is also realized in 
the simple (ffc) hopping model on an fee lattice even 
without oxygen p orbitals^ These findings indicate that 
the tendency toward antiferro-octupole ordering is 
common to Tg models on fee lattices. 
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